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We propose a constituent quark model to evaluate heavy decay constants and form factors relevant 
for B D^*^ semileptonic transitions. We show that the model reproduces the scaling laws dictated 
by the spin-flavour symmetry in the heavy quark limit and describes quite well the experimental 
data. 



I. INTRODUCTION 

■ n 

' The study of exclusive charmed semileptonic decays of B mesons is of primary importance to extract |1| one of the 
. free parameters of the Standard Model: the absolute value of the Cabibbo-Kobayashi-Maskawa matrix element, | Vcb| 
The extraction is based on the prediction of the Heavy Quark Effective Theory (HQET) j3| which fix an absolute 
, normalization, at zero recoil point, of the form factor which survives in the limit of infinite quark masses. Moreover, 
5-) ' it is possible to show that differently from B — > Diiy process the B — > D*£i/ decay doesn't receive l/mg corrections 
' at zero recoil point 0|. This facts allow us to extract \Vcb\ from the differential partial decay width for the B —> D*£iy 

^ process in a nearly model independent way p, . 

In this paper we will study the B D^*'ii' processes from a different point of view. We propose a very simple 
constituent quark model to evaluate heavy decay constants and heavy-to- heavy form factors. They exhibit the scaling 
T-H ■ laws dictated by the HQET at leading order and describe in a satisfactory way the experimental data. To study 
^ , the semileptonic transitions between the heavy mesons B and Z?^*-' and to compute the relevant hadronic matrix 
• elements we use the ideas presented in the papers in Ref. 7] devoted to study heavy to light semileptonic and rare B 
] transitions. In that papers, the heavy meson B is described as a. bq {q £ {u, d}) bound state and the corresponding 
' wave functioii, '(/'(fc), is obtained by solving a QCD relativistic potential model. Here, we adopt a different point of 
I view. As in Q, we describe the involved (heavy) mesons as bound state of a heavy quark and a light anti-quark 
I but for the wave functions we assume their mathematical form and we fix the free parameters by comparing model 
predictions and experimental data, when available (see after). 

The paper is organized as follows. In the next section we introduce our constituent quark model, heavy decay constants 
O • and heavy-to-heavy form factors are evaluated in section IIIII The section IIVI is devoted to discuss the heavy quark 
I " limit for decay constants and form factors. Numerical results and conclusions are collected in section Ivl 

D : 



II. THE MODEL 



We describe any heavy meson H{Qq), with Q G {6, c}, by introducing the matrix 



- 1 — \ 

X 



H^-r=i^H{k)j ^--^ — ^r— ^, (1) 

V3 V ^q''^Q +91-92 2mQ 2mq 

where toq (m^) stands for the heavy (light) quark mass; g^, q!^ their corresponding 4— momenta (cf. Fig. With 
ipnik) we indicate the meson's wave function and the factors are chosen to satisfy the following relations 

< H\H > = Tr{(-7oi?t7o) H} = 2 mn , 

^|V./(fc)P = 2mH. (2) 

The meson-constituent quarks vertex, F, is given by 



F = — «75 = Fp for pseudoscalar mesons (3) 

E Fy(£, gi, (72) for vector mesons (4) 



rriH + mq + 



where e is the polarization 4— vector of the (vector) meson H . In any HQq vertex we assume the 4— momentum 
conservation, i.e. 9i +92 = P^^ ^he H meson 4— momentum. Therefore, if we choose — {Eg^k), and q!^ — (Eg, — fc). 
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i.e. the H rest frame, we have {k = \k\) 

Eq + E^ = ^m^ + k^ + ^ml + k^ = niH ■ (5) 

which can be read as the definition of a running heavy quark mass, as was done in Ref. 0. In fact, the Eq. (jsj with 
the constraint mQ(fc) > gives the relation 

? 9 

0<k<KM = — f (6) 



on the loop momentum k 

d^k 



(7) 



Let us now write down the remaining rules for the computation of the hadronic matrix elements in the framework of 
this model: 

a) for the weak hadronic current, ^2 T'^ qi, one puts the factor 



where is some combination of Dirac matrices; 

b) for each quark loop, in addition to the integration in Eq. |(7J), one puts a colour factor of 3 and performs a trace 
over Dirac matrices. 

III. LEPTONIC DECAY CONSTANTS AND B -> D^*^ SEMILEPTONIC TRANSITIONS 

In this section we introduce heavy decay constants and semileptonic form factor for heavy-to-heavy transitions and 
we give their expressions in the framework of our model. 

Using the rules introduced in the previous section we immediately get the expressions for the heavy meson decay 
constants. The pseudoscalar case was obtained and discussed in Ref. 0, for future convenience, we report the 
resulting expression for the B meson: 



V3 f ,,,2, ,,Amb + mq){mbmq+qi-q2) 
fB = o 2 2 / dkk il)B{k) ^ ■ (9) 

Moreover, we have evaluate the vector heavy meson decay constant, which is defined by 

< Q\V^\H*{p,e) > = TOff. /h. e^. (10) 
In particular, if we consider the B* meson, we obtain 



V3 



/s* = TT^ / dk 



271^1713' Jo ^ EbEqinibing + qi ■ (72) 



{mbniq + qi- q^) - - 



2 k^ruB 



3 viB' + mb + m„ 



(11) 



B ^ D and B — > D* form factors 



The same rules allow us to evaluate the matrix element < D{p')\cjp,b\B{p) > relevant to the weak semileptonic 
transition of B to D mesons. With reference to the graph in Fig. ^and choosing the 4— momenta qi and 92 as in the 
previous section and qt^ = (Ec, k — q), we get 



< D{p')\cj^b\B{p) >-- 



Sk 
(2^ 



ipD{k)ipB{k) 



mq-mb 



rnqTric 



EbEc 



Tr 



'2 +™9/-pt +nT-c il 



2m, 



2m, 



rnqTUb + gi • 92 V "'-'j'Tlc + 93 • 92 

"^b /T- A ^ ^2 + rriq 
(1 p) 



2m„ 



(12) 
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FIG. 1; Quark model diagram for the semileptonic B decays involving b ^ c transition. The thin lines represent quarks, the 
thick ones mesons. The gray disks represent the quark-quark-meson vertexes. 



In the previous equation the integration domain V is fixed enforcing the energy conservation both in the initial and 
final quarks-meson vertexes. This can be done introducing, in addition to the beauty running mass (cf. Eq. l(SJ)), the 
charm running mass mc(fc) for which mc(fc) > 0. After some algebra the physical domain V is found to be given by 



Max(0,A:_)< k <Mm{KM,k+) 
Max(-l,/(fc,|g|)) < cos{9) < +1 
< (h < 27r 



with 



fik, k1) 



|g1 {mjj + ml) ± {m% ~ m^)^TO|, + q' 
2m|, 

2^^ ■ 



(13) 



(14) 



(15) 



(j) and 9 are the azimuthal and the polar angles, respectively. Note that we have chosen the z— axis along the direction 
of q, the (tri-)momentum of the W boson (cf. Fig. 

The Eq. (|12|l allows us to immediately extract the form factors f±{q^) defined by 



(16) 



The last matrix element relevant to charmed semileptonic decay of B mesons is usually written in terms of the 
following form factors 

<D*ip',s)\cj^il-j5)b\Bip) > ^2g{q') e^^^^^ elp^p'f, 

[ a+('z') (Pm+p;) + (Pm-p;)]} , (17) 

they are connected in our model to 

d^k 



<i5*(p',£)|c7p(l"75)fo|i?(p) > = 



fo \i'^D'{k)ipB{k). 



EhEr 



Tr 



2m„ 



-(rv(e,g3,g2)t)^;^7.(l-75)-2^^^ 



2to„ 



(18) 



Also in this case the extraction of the form factors in Eq. ()17|l can be done using the same frame we adopt for the 
extraction of /±(q^). For the polarization vectors we use: 



(0,-1,0,0) A = l 

e^{\)^{ (0,0,1,0) A = 2 

(|g|,0,0,-£;D.)Mc* A = L 



(19) 



where Etj- (= ■^/(f ^ + mjjt) represents the energy of the D* meson. 
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IV. HEAVY QUARK LIMIT 

In this section we discuss the Heavy Quark Limit for decay constants and form factors. We show that decay 
constants and heavy-to-heavy form factors satisfy the scahng laws predicted by HQET at leading order jj] . 
To show how the results of our model depend on the heavy quark mass, we need to specify the shape of the wave 
functions tpnik). We choose two possible form, the gaussian-type, extensively used in literature (see for example 

Si) 



, ,3 ^^i' 1 9 , ,2 
H L H 



i^H (k) = 47r3/4 /!!^ exp <i , (20) 



V'//(fc) = 47r./^exp(— ^ (21) 



and the exponential one 

which is able to fit the results of relativistic quark model Td]. In our approach ujh is a free parameter which should 
be fixed (cf. next section for details). 

A. Heavy Decay Constants 

To extract the heavy mass dependence from the decay constant, it is useful to define x = {2ak)/mB in such a way 
the expressions in Eq. lO and (|ll|l can be formally written as 

/•a 

fBM= dx yjB{k{x)) Fbm{x,z), (22) 

JQ 

where Fb(-) (x, z) have a very simple expressions for 2 = (z = mq/niB) 

V 2 Stt^ a-^ ^[a - x){2a - x) 

™2 „2 3 (a + ^y a(a — x)] — x 

Fb4^,0) - } \ r- ?^=T (24) 

8v6 TT-^ a-^ \'2a — X [y/a + ^/a ~ x) 



The integral in Eq. H22|) , for < a ^ 1 can be evaluated analytically, obtaining for the leading behaviour the following 
result 



1 \/6uj% 

gaussian — type 



" ^ 1 . , (''^ 

exponential — type 



in both cases in agreement with the scaling law predicted by the HQET. 

B. B ^ D Form Factors 

The same procedure applied to the heavy-to-heavy {B D) transitions allows us to find the scaling laws of the 
form factors f± defined in Eq. (|16|) . As for decay constants, we can formally write 



f±{qn = / dx ^B{k{x)) ^JD{k{x)) F±{x, z, q') , (26) 
Jo 

where, for z = 0, a; ^ 1 and near the zero recoil point (q^ = q^y^^^^) 

F±{x,0, g')|,^^,^„^ ^ -^^^ mlimD ± tub) U - ^^{w - l)\ . (27) 
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Here w = v -v' with v and v' the four-velocities of the B and D mesons, respectively. Also in this case we can extract 
the dependence of the form factors from the heavy masses performing the integration in Eq. (|2(j|l 



LOB 



cut 



2^m.DmB 



2V2 

{lob + ^dY 



LjJt 



3/2 



11 

12 



11/ 

i-Y^(--i) 



gaussian — type 



exponential — type 



(28) 



It should be observed that the terms in square brackets should be interpreted as the Isgur-Wise function, ^(w), near 
to ti; = 1. Moreover, in the Heavy Quark Limit we should have lob — (^d which implies the correct normalization, 
^(1) = 1, for both wave functions. 



C. B ^ D* Form Factors 



The same analysis can be carried out for the B D* form factors. Let us start to consider the form factors a±{q^) 
As for the f±, we can write 



dx ipBiHx)) ipD^mx)) A±{x, z, q^) , 



where, for z = 0, a: ^ 1 and near the zero recoil point 

A±{x,Q,q^)\q2^q2^^^ 



--(l-i(^-l) 



(29) 



(30) 



QA'K'^Ot' TUB 

Analogously to the B ^ D case, the heavy mass dependence can be obtained performing the integration in Eq. I|29|l . 

3/2 



a±(9^)lg2~g2 



1 



lol+Lol 



2^mDmB 



2V2 



11, 

i-Y^(--i) 



[UJB + OJD ) 



11 / 

1 iw - 1 

12^ 



gaussian — type 



exponential — type 



(31) 



The behaviour with heavy masses of the vectorial form factor, g{q^), is the same of a^{q'^) in agreement with the 
prediction of Heavy Quark Symmetry. For the last axial form factor, f{q^), our model predicts 



/(9^)lg2~g2^^^ y/niD-mB (1 + w) < 



2V2 



OJB UJD 



3/2 



11 . 

l-^(--l) 



1- — (w-1) 
12^ ^ 



gaussian — type 



exponential — type 



(32) 



Thus all the form factors satisfy the scaling laws dictated by the HQET. Moreover, the model predicts the following 
Isgur-Wisc function: 



i(w) = 1 - ii(^ - 1) + g(^ - 1)2 + o{{w - If) 



(33) 



where the quadratic term, neglected in Eqs. (|28|) . H31|l and (|32|l . is shown. The resulting Isgur-Wise function satisfies 
both the Bjorken Sum Rule hj 



P^--e'(l) = ^>7. 



11 ^ 3 
12-4 



(34) 



and the lower bound on the curvature Il2 



-^-r(i) = 5>^P^(i + >^) = 



99 

80 



(35) 
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Exp. or Lattice 


our fit (exp.) 


our fit (gauss.) 


fD/fB 


1.23 ±0.22 [13] 


1.03 


1.05 


Id 


300ll«°l«o MeV fl] 


145 MeV 


145 MeV 


BR(B Dlu) 


(2.14 ±0.15)% fl] 


2.04 % 


1.75 % 


fB 




140 MeV 


139 MeV 



TABLE 1: The experimental values |1||] and Lattice result for the decay constants used in the fit of the free parameters of 
the model. For the free parameter we assume ujd = <^d*. Moreover, we use \Vcb\ = (41.3 ± 1.5) X 10"^ Q. 



Parameter 


fitted values (exp.) 


fitted values (gauss.) 


ruq 


311 MeV 


269 MeV 


UlB 


258 MeV 


421 MeV 




255 MeV 


347 MeV 



TABLE II: The values for the free parameters of the model in correspondence of the best fit for both wave functions. The two 
sets of values are obtained for the exponential (exp.) and gaussian vertex (gauss.). 

V. NUMERICAL RESULTS AND DISCUSSION 

As we have seen in section Hill the heavy-to-heavy form factors can be easily extracted with the help of the Eqs. H12|l . 
(|15|l . H17|) and (|18|l . Nevertheless, unlike for the decay constants, their analytical expressions are quite long, and, for 
the sake of brevity, we do not report them here. 

As already discussed in the previous section, to evaluate numerically form factors and decay constants, we must 
fix the meson wave functions, ipuik). For the wave function we considered two possibilities: the gaussian and the 
exponential form. In both cases, for any heavy meson, H, we have one more free parameter, luh- In order to determine 
the free parameters of the model we proceed as follows. We neglect differences between pseudoscalar and vector mesons 
in the vertex function, in other words we put lod = wd* • Moreover, we neglect differences between u and d quark 
masses. In such a way the free parameters of the model are lub, t^D and rriq. They are adjusted by fitting the 
experimental values of fu, BR(_B — s- Div) and the results of lattice simulation on the ratio fo/fs- The numerical 
results are collected in Tables III and Ull 

Comments about the results in Table ^ are in order. Let us start with decay constants. The model predicts 
large l/rric corrections for in such a way the predicted ratio fo/ fs violates strongly the heavy quark mass limit. 
However, in the allowed region for the parameters there is the possibility to fulfill both the heavy quark limit and the 
Lattice result but the values of the decay constants [fo ^ fs ^ 140 MeV) are predicted smaller than the ones obtained 
by Lattice simulations Regarding the B ^ D form factors, it should be observed that the experimental value for 
the BR(_B Dlv) can be reproduced with \Vcb\ ~ (41.3 ± 1.5) x 10~^ j^]- However, when the exponential function is 
considered, the agreement becomes better as a consequence of the larger value predicted for /+(0) (/+(0) = 0.57 (0.51) 
for exponential (gaussian) vertex function). The same situation occurs if the differential partial decay width for the 
B D*tv is considered. In FigureEl assuming the values in Table ITTI we plot dT{B D*£i/)/dw in comparison with 
experimental data ^,1^. In particular, the left (right) panel allows to compare model predictions and experimental 
data for the exponential (gaussian) vertex function. Both panels contain three curves corresponding to the predicted 
dT(B D*(iv)/dw for the central, upper and lower l-cr values of \Vcb\ Q- The agreement between model predictions 
and experimental data is quite good, a better agreement requires a smaller value of \Vcb\- In this respect, using 
exponential vertex function and a value of \Vcb\ = 38.7 x 10~^ 0, the predicted dT{B D*£h')/dw is in very good 
agreement with experimental data from Babar Q. 

In conclusion we have proposed a constituent quark model to describe heavy mesons. We showed that the model 
predictions on decay constants and form factors reproduce the scaling laws dictated by HQET at leading order in the 
heavy quark mass limit. For finite heavy quark masses the agreement with experimental data is quite good. 
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FIG. 2: Predicted ranges for dT{B D*£v)/dw compared to data. Solid boxes (triangles) refer to B° ^ D+*rv {B' ^ 
D°*£~i') process l5|. Data points from BABAR are displayed with stars. Solid lines refer to model predictions for exponential 
(left) and gaussian (right) in correspondence of |Vc6| = (39.8,41.3,42.8) x 10~^ Q- 
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